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Abstract. Dispersive Alfv6n wave activity is identified in 
four dayside auroral oval events measured by the Freja satel- 
lite. The events are characterised by ion injection, bursty 
electron precipitation below about 1 keV, transverse ion heat- 
ing and broadband extremely low frequency (ELF) emissions 
below the lower hybrid cutoff frequency (a few kHz). The 
broadband emissions are observed to become more electro- 
static towards higher frequencies. Large-scale density de- 
pletions/cavities, as determined by the Langmuir probe mea- 
surements, and strong electrostatic emissions are often ob- 
served simultaneously. A correlation study has been carried 
out between the E- and B-field fluctuations below 64 Hz (the 
dc instrument’s upper threshold) and the characteristics of the 
precipitating electrons. This study revealed that the energisa- 
tion of electrons is indeed related to the broadband ELF emis- 
sions and that the electrostatic component plays a predomi- 
nant role during very active magnetospheric conditions. Fur- 
thermore, the effect of the ELF electromagnetic emissions 
on the larger scale field-aligned current systems has been in- 
vestigated, and it is found that such an effect cannot be de- 
tected. Instead, the AIfv6nic activity creates a local region of 
field-aligned currents. It is suggested that dispersive Alfv6n 
waves set up these local field-aligned current regions and in 
turn trigger more electrostatic emissions during certain con- 
ditions. In these regions ions are transversely heated, and 
large-scale density depletions/cavities may be created during 
especially active periods. 


1 Introduction 

An AlfV6n wave propagating from the magnetosphere to- 
wards the ionosphere will encounter a denser plasma and 
as a result may change its characteristics. For instance, the 
Alfvgn wave may be reflected (e.g. Haerendel, 1983; Knud- 
sen et al., 1990, 1992), dissipated (e.g. Gark and Seyler, 
1999; Goertz, 1984), or broken down into other wave modes 
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(e.g. Seyler et al., 1998). A dispersive Alfv6n wave (DAW) 
(e.g. Lysak and Lotko, 1996; Stasiewicz et al., 2000, and ref- 
erences therein) may therefore enhance both its inertial and 
kinetic characteristics, and develop a field-aligned electric 
field in the wave front when it propagates at a small angle to 
the Earth’s magnetic field. These AIfv6n waves are therefore 
often suggested as candidates for auroral zone field-aligned 
electron acceleration process. 

The low frequency part of the broadband extremely low 
frequency (ELF) emissions observed by spacecraft in the au- 
roral zone is often considered to be Alfv6nic in nature, with 
perpendicular wavelengths of the order of the electron iner- 
tial length (A e ) (e.g. Chmyrev et al., 1989; Louam et al., 
1994; Wahlund et al., 1994; Seyler et al., 1995), and to be 
associated with field-aligned electron bursts at suprathermal 
energies (Boehm et al., 1990; Knudsen and Wahlund, 1998; 
Wahlund et al., 1998; Chaston et al., 1999). This wave activ- 
ity is also connected with transverse heating of ions (Knud- 
sen and Wahlund, 1998; Knudsen et al., 1998) and might 
cause large-scale density depletions. 

The field-aligned currents into and out from the ionosphere 
are often described as being closed within the ionosphere. 
The conductivity in the ionosphere is in this scenario con- 
nected to the strength of the field-aligned current and hence 
the particle acceleration. When a widespread region of DAW 
activity exists, this scenario will change to a local region of 
very strongly fluctuating field-aligned currents of several 100 
pAJm 2 with a rather weak contribution to the overall large- 
scale current systems. 

The Freja mission is dedicated to observe the space plasma 
environment in polar regions with an orbit altitude of 900- 
1700 km. In this paper we apply the formalism of the con- 
cept of DAW (e.g. Lysak and Lotko, 1996) to the ELF electric 
and magnetic field emissions measured by the Freja space- 
craft in order to show that these measurements are consistent 
with such a model. Other observed parameters, such as those 
associated with field-aligned electron bursts, transverse ion 
heating, field-aligned current systems and the appearance of 
large-scale plasma density depletions/cavities, are discussed 
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in terms of a DAW model. In sections one and two the char- 
acteristics of four Freja passages through the dayside auroral 
region are presented and discussed. In section three a de- 
tailed study of the ELF emissions below 64 Hz is presented 
and discussed in terms of a DAW model. In section four 
the expected field-aligned electric fields and propagation an- 
gles are derived from the electric and magnetic field measure- 
ments under the assumption that they are produce by disper- 
sive Alfv6n waves. In addition, the DAW associated field- 
aligned currents are investigated with regard to the effect on 
field-aligned current systems. A correlation study between 
the electric and magnetic ELF fields and field-aligned elec- 
tron burst characteristics is described, and finally we make 
some concluding remarks regarding the overall interpretation 
of the observed processes in terms of wave-particle interac- 
tions caused by dispersive Alfv6n waves. 

2 The Freja instrumentation and orbit characteristics 

The Swedish/German satellite Freja passes the northern au- 
roral region at an altitude of about 1700 km and at an inclina- 
tion of 63°. It is a sun-pointing spacecraft with a spin period 
of 6 seconds. Detailed information about the Freja experi- 
ments can be found in a special Freja issue of Space Science 
Reviews (see Boehm et al., 1994; EJiasson et al., 1994; Hol- 
back et al., 1994; Marklund et al., 1994). Relevant experi- 
mental information for this study is discussed where neces- 
sary. 

3 General characteristics of the four Freja auroral oval 
crossings 

Data are presented from four Freja satellite passages through 
the dayside auroral region during ion injection events. For 
convenience these orbits are referred to as /, 77, III , and 
IV , respectively, corresponding to actual orbits 5200 (Novem- 
ber 3, 1993), 5265 (November 8, 1993), 5292 (November 10, 
1993), and 6653 (February 21, 1994). A short description 
of the general observational characteristics is given here, and 
more detailed information regarding these particular data can 
be found in (Andersson et al., 2000). 

Figure 1 displays the overview of particle and wave data 
from all four orbits, and Figure 2 shows the corresponding 
energy densities as well as the dc electric and magnetic field 
waveforms. The three first orbits occurred during moderately 
active conditions (K p ~2 — 3), while orbit IV occurred during 
very active conditions (K p ~6). 

3. 1 Ion injection and heating 

All four orbits are characterised by clear ion injections with 
energies from 500 e V up to the TICS instrument threshold of 
4.3 keV (Figure 1, panel b). The downward proton energy 
flux (Figure 2, panel a) is calculated by integrating over all 
TICS energies from 50 eV to 4.3 keV and over pitch angles 
between 0°and 90°. The time resolution of this calculated 


ion energy flux is 400 ms. Occasionally these fluxes are 
somewhat spin modulated due to limitations in pitch-angle 
coverage. Downward ion energy fluxes of 0.01-1 mW/m 2 
exist during the observed ion injections. The ion injections 
are accompanied by bursty electron precipitation, transverse 
ion heating as well as broadband ELF emissions up to the 
lower hybrid cutoff frequency (Figure 1). Large-scale den- 
sity depletions, as measured by the Langmuir probe, are only 
detected during the most active orbit, TV (Figure 1 , panel 6). 
Intense transverse ion heating (mainly 0 + ) is only encoun- 
tered in orbit IV, for which it reaches energies of up to ~400 
eV, while the ion heating during the other three orbits reaches 
at most ~50 eV. 

3.2 Electron precipitation characteristics 

The electron bursts have downward energy fluxes (calculated 
from pitch angles 0°-90°) of the order of 0.01-10 mW/m 2 
(Figure 2, panels b and c), with the most intense events en- 
countered in orbit IV. Panel b contains TESP data in the 
lower energy range, 12-50 eV, which corresponds roughly to 
energies below the local Alfv6n velocity (tu). Panel c in- 
cludes almost the full energy range of the TESP instrument 
(36 eV-25 keV). The time resolution of the electron mea- 
surements is 62.5 ms. 

The most field-aligned electron precipitation can also be 
found in orbit IV (Figure 2, panel d; red line denotes 0° pitch 
angle, green 45°and black 90°). In all four orbits the flux 
in 45°and 90°pitch angle is equal but during the bursts of 
electrons the 0°pitch angle fluxes are much higher. In orbits 
II and III bursts appear more sporadically, while in orbit 
IV the electron bursts occur more frequently in time. 

3.3 Cold plasma characteristics 

A set of spherical Langmuir probes samples continuously the 
current from the plasma at a fixed positive bias voltage. This 
current is proportional to n c /\/T7, where n t and T e are the 
electron density and temperature. Usually the variation of 
the electron temperature is small and slow and the Langmuir 
probe current variation is assumed to be due solely to plasma 
density variations. When possible, we have used the nar- 
rowband HF Langmuir emissions near the plasma frequency, 
u pe , as well as the intermittent Langmuir probe bias voltage 
sweeps to confirm that the sampled probe current is mainly 
due to the variations in n e , even though a variation in T c may 
account for a current decrease by up to a factor of 2 during 
active periods. 

The regularly measured probe sweeps are seen as shaded 
grey areas in Figure 2 (panel e). In Table 1 measurements 
from the analysed probe sweeps, such as the plasma temper- 
atures (Ti, T e ) and the temperature ratios (T,/T c ), are pre- 
sented together with some useful calculated plasma charac- 
teristics. As mentioned above, the ion heating is weak during 
orbits II and III , somewhat stronger in orbit I and intense 
in orbit IV. It should be noted that it is the whole bulk ion 
population that is heated, not just the energetic tail of the ion 
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Fig. 1. Overview data from four Freja satellite passes of the northern dayside auroral region. The panels display from the top; 0+ (TICS), H+ (TICS), 
field-aligned downgoing electrons (TESP), HF electric wave component (F4), MF electric and magnetic wave components (F4), and the ac magnetic field (F2) 
respectively. 






L. Andersson et al.: Electromagnetic waves and bursty electron acceleration 



TIME (seconds) TIME (seconds) 

Fig, 2. Data from the same four Freja satellite passes as shown in Figure 1. Hie panels show from the top; proton (50 eV-4.3 keV), electron (12-50 eV) and 
electron (36 eV-25 keV) downward energy fluxes (mW/m 2 ); the total electron counts for dace different pitch mgles (0°rcd, 45° green nd 90° black); the 
Langmuir probe measured density; the dc electric and magnetic fields perpendicular to the model magnetic field lines. 
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distribution. 

During orbit TV, several density drops (Figure 2, panel 
e) are observed simultaneously with ion injection (panel a) 
and high electron count rates in the field-aligned direction as 
compared to count rates at other pitch angles (panel d). The 
three largest density drops (hereafter referred to as cavities) 
are centred around +55, +90 and +135 seconds into the plot. 
The highest electron count rates ( 1 20- 130 seconds) have their 
main contribution from the lowest energies, below 50 eV. The 
electric field fluctuations (panels f and g) during orbit IV 
have extreme amplitudes (>3 V/m!), especially within the 
cavity regions. 

3.4 Electromagnetic fields and wave characteristics 

It is clearly seen from Figure 1 (panels e and f) that the broad- 
band ELF emissions can reach frequencies above 1 kHz and 
in some cases up to 6 kHz, which is close to the lower hy- 
brid cutoff. It can also be seen that the emissions become 
increasingly electrostatic towards higher frequencies, while 
the magnetic component becomes dominant at the lowest fre- 
quencies (see panels f and g, for ac magnetic field measure- 
ments). See Wahlund et al. (1998) for further details regard- 
ing the broadband ELF emissions. This study concentrates 
more on the lower frequency, more electromagnetic part, of 

ELF spectrum and its relationship to electron bursts, day- 
side ion injections and density depletions. 

The dc electric and magnetic field instruments were op- 
erated somewhat differently during the four selected orbits. 
The dc electric field measurements normally make use of two 
orthogonally positioned spherical probe pairs on 10.5 m long 
wire booms in the spin plane (i.e. 21 m between spheres). 
However, during orbit IV one of the probe pairs was oper- 
ated in a current collection mode allowing for 6n/n interfer- 
ometry measurements, and as a consequence only one probe 
pair measures the dc electric field on this orbit. This explains 
why the electric field on orbit IV (Figure 2, panels f and g) is 
spin modulated. The sampling rates for orbits 7, II and III 
are 768 samples/s for the dc electric field, and 128 samples/s 
for the dc magnetic field. On orbit IV the dc electric field 
was sampled at 1536 samples/s and the dc magnetic field at 
256 samples/s. 

The dc electric and magnetic field measurements have been 
de-spun and are presented in a coordinate system for which 
the x-axis is along a model magnetic field line, positive y 
points towards west, and positive z points towards the equa- 
tor. The Freja satellite measures the electric field in the spin 
plane of the spacecraft and therefore only gives a 2-dimen- 
sional estimate of the total electric field. We therefore as- 
sume that the magnetic field-aligned component (E x ) is zero. 
The assumption should be good since the field-aligned elec- 
tric field component of a dispersive Alfv6n wave is theoreti- 
cally much smaller than the perpendicular electric field com- 
ponent. 

The field fluctuations (Figure 2, panels e to f) are of ex- 
tremely large amplitude during orbit IV. When the down- 
ward ion energy flux is close to 1 mW/m 2 the magnetic field 


data show a large gradient, the electron distribution is isotropic 
with a loss cone, there is no density gradient, and the dc elec- 
tric field shows just a spin modulation about zero. On orbit 
III downward bursts of electrons are correlated with sharp 
changes of the magnetic field. 


4 Detailed study of the ELF emissions below 64 Hz 

The time series of the dc electric and magnetic field mea- 
surements were analysed by an overlapping sliding window 
Fast Fourier Transform (FFT) analysis (Figure 3, panels a-e). 
In order to cover all frequencies below 64 Hz, four different 
time windows were selected (15, 6, 1, and 0.3 s correspond- 
ing to lowest non-dc frequencies of 0.133, 0.333, 2 and 7.1 
Hz respectively). Each time window was weighted with a 
Hamming filter. The longer time windows are then used at 
the lower frequencies, while the shorter time windows are 
used at the higher frequencies in Figure 3 (panels a-e). The 
wave spectrograms presented in Figure 3 have the same ap- 
parent time resolution as the electron data (TESP, 62.5 ms). 
This analysis corre sponds very much to a wavelet analysis, 
although it is FFT-based. 

4.1 Poyntingflux 

The wave energy per unit time and area is represented by the 
spectral Poynting flux (S (/) = \/^B (/) x B(/)) in Fig- 
ure 3 (panel f). The Poynting flux along the Earth’s magnetic 
field is calculated from the wave data represented in panels a 
through e by 

S * = te[j-(E V B Z -E Z B„)} (1) 

for each frequency component. It should be noted that the 
direction and magnitude of the total Poynting flux when inte- 
grated over frequency will depend on the frequency window 
used and their significance is unclear when treating broad- 
band plasma wave emissions with several possible non-linear 
wave modes involved. We do not present the direction of the 
spectral Poynting flux because it shows very complex small 
scale variations, more suitable for detailed event studies on 
the time scale of seconds, and is therefore left for a future 
paper. Here we just note that the spectral Poynting flux is 
highest at the lowest frequencies, and correlates well with 
the downward electron energy flux (Figure 2, panels b-d). 

4.2 The Alfv6n velocity and the 5E± /5B±-ra tio 

Panels g and h in Figure 3 display as function of time the ra- 
tio between the electric and magnetic field fluctuations 
(SE±/6B±) and the local Alfv6n velocity (t^). The local 
Alfv6n velocity == Bo/ y/iMUTiiTii) is calculated for four 

different ion compositions (using a H + and 0 + mixture with 
100%, 88%, 47%, 0% composition) and is based on the 
plasma density estimated from the Langmuir probe measure- 
ments. 
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Fig. 3. More detailed ELF wave characteristics for the four orbits presented in Figure I . Panels from the top are respectively; a-e) the power spectral densities 
(PSDs) for the three magnetic field components and the two perpendicular electric field components; f) the magnitude of the Poynting flux (5); g) 5E±/5Bj_; 
h) the AIfv£n speed (va) for four different plasma compositions; i-j) the value of A e k± for a pure oxygen and a pure hydrogen plasma respectively; k) the 
parallel wave vector (fey); I) the parallel group velocity of inertial Alfv6n waves as fraction of the local Alfvdn speed; m) the perpendicular group velocity of 
inertial Alfv6n waves. See the text for more details. 
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Table L Values fr om the probe sweeps 


ORBIT 

/ 

II 

II 

II 

III 

III 

IV 

IV 

N 0 +/N ittot 

0.85 

0.95 

0.95 

1.0 

0.8 

0.9 

0.9 

0.95 

T t [eV] 

0.8 

0.7 

0.7 

0.6 

0.5 

0.6 

0.5 

2.5 

T, [eV] 

0.5 

0.1 

0.1 

0.2 

0.1 

0.1 

0.1 

15 

N t [cm" 3 ] 

6100 

5000 

5000 

4200 

3000 

2500 

5000 

4000 

Ti/Te 

0.63 

0.14 

0.14 

0.33 

0.33 

0.17 

0.2 

6 

Bo [nT] 

29400 

29900 

29700 

29500 

29500 

29400 

30600 

30800 

Po+ M 

9.8 

4.3 

4.4 

6.2 

4.4 

4.4 

4.2 

51.4 

[m] 

68 

75 

75 

82 

97 

106 

75 

84 

va [10 6 m/s] 

2.1 

2.3 

2.3 

2.5 

2.9 

3.2 

2.4 

2.7 

CGLAT [°] 

74.5 

73.4 

73.8 

74.2 

73.5 


67.3 

70.5 

H [km] 

1730 

1660 

1700 

1720 

1710 

1725 

1537 

1570 

UT 

18:47:10 

16:45:05 

16:47:25 

16:49:45 

17:49:00 

17:51:50 

17:01:30 

17:03:40 


Data in the top five rows originate from evaluation of the probe sweeps. The values in the middle four rows are calculated using the assumption of a 100% 
oxygen plasma. The last three rows show the locations and times of the spacecraft during the probe sweeps. 


At the lowest displayed frequencies the ratio 6Ex/SBx 
is lower than va during the three first moderately active or- 
bits (J, II and III). This suggests that dc-like field-aligned 
current systems dominate the magnetic field fluctuations be- 
low about 0.5 Hz. Between 0.5-7 Hz the 5E±/6Bx -ratio 
increases to a more stable plateau value somewhat above the 
local Alfv6n velocity for a pure oxygen plasma. The drop at 
the highest displayed frequencies is probably an instrumen- 
tal artifact; it might be that the magnetic emissions reach the 
noise level of the fluxgate magnetometer. 

During orbit IV the whole displayed frequency interval at- 
tains SEx / & B ± -ratios in excess of the local Alfv6n velocity 
for all possible H + -0 + compositions. The ELF emissions 
are therefore more electrostatic than expected for a classical 
Alfv6n wave, and inertial and/or kinetic effects are probably 
important. The large values at the lowest frequencies below 

0. 5 Hz may just be an effect of that only one probe pair was 
used for electric field measurements. 

4.3 Dispersive Alfvdn Waves (DAW) 

The term DAW is a general description covering inertial 
Alfvdn waves (I AW, Goertz an d Boswell (1979)) applica- 
ble in a plasma with v t h, e (= y/kBT e /m e ) < va and /?(= 
2pon e ksT e /B^) < m c /mj, and kinetic Alfv6n waves (KAW, 
Stefant (1970); Hasegawa (1976)), applicable in a plasma 
with t kh,e > va and /3 > m e /mi. The concept originates 
from a development by Lysak and Lotko (1996), and is de- 
scribed in more detail in a review article by Stasiewicz et al. 
(2000). The dispersive properties of DAWs become impor- 
tant when the perpendicular wavelength becomes compara- 
ble to the ion gyroradius (pi) or the electron skin depth (A c ), 

1. e. when the quantities kxPi and/or kx A e become compa- 
rable to or larger than 1. When this happens parallel electric 
fields develop, which in turn can accelerate electrons. 

The dispersion relation for a low frequency DAW (u < 
u> C i) propagating in a cold electron (t^e < va) low (3 plasma 


can be described 

2 _ jl2 2 (1 + 

" -v* 0 + 4 *?) 

and the group velocity becomes 


dw 


VA 


-ywA e 


k±Xe 
1 + Jfc’A? 


( 2 ) 


(3) 


dk X (1 + *3A2)i/2 

Solving the dispersion relation in this cold electron limit 
gives a ratio of 


I ^ 1= V A \J{\ + A= ) ( 1 + 


(4) 


and the relationship between the field-aligned and the per- 
pendicular electric field becomes 


El = k jM m 

E ± 1 + *2 XI W 

These relationships do not take into account several possible 
non-linear effects that can for example, steepen the wave- 
front of a DAW, breaking the wave into higher frequency 
components, and enhancing the parallel electric field further 
(e.g, Seyleret al., 1995; Seyler and Wahlund, 1996). 


4.4 Inertial term and kx 


The perpendicular wave vector (kx) can be estimated from 
the Alfv6n speed (r^) and the ratio SEx/SBx (equation 4). 
From Table 1, p? << A^, which is a good approximation 
except possibly during the active periods on orbit TV. The 
value of kx A e is presented in Figure 3 (panels i and j) for 
pure hydrogen and oxygen plasmas respectively. It should 
be noted that only points which satisfy SEx/SBx > va for 
a pure oxygen plasma are kept in panels i through m, and 
values below va are white in those graphs (equation 4). 

The value of kx A e increases significantly above zero when 
enhanced broadband electric field emissions appear (panels d 
and e), and values up to kx A e ~l-2 can be detected during 
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orbit IV. This correspond to perpendicular wavelengths of 
around 250-500 m in the case of orbit IV and values in ex- 
cess of that for the other three orbits. 

The maximum possible Doppler shift (broadening) of the 
perpendicular wave vector (wave fields) due to the motion of 
the satellite is estimated to be k±v tc ~ 10-20 Hz for orbits /, 
II and III and 100-200 Hz for orbit IV. Thus the enhanced 
ELF wave emissions which occur during orbit IV may be 
heavily Doppler broadened (shifted) depending on the exact 
k-vector distribution, 

4.5 Group velocity and k\\ 

From knowledge of k ± , several other characteristics can be 
estimated; among them the perpendicular and the parallel 
group velocities (panels 1 and m) as well as the parallel wave 
vector (panel k). These parameters are estimated assuming a 
pure oxygen plasma and by using k± calculated from panel j . 
While k± is independent of frequency, except for frequencies 
above about 40 Hz during active parts of orbit TV , will de- 

pend on frequency throughout the DAW dispersion relation. 
The parallel wavelengths as estimated in this way are in the 
range of 100-1000 km. 

The parallel group velocity is close to the Alfv6n speed 
during orbits J, II and III , while during orbit IV the par- 
allel group velocity becomes a fraction of the Alfv6n speed, 
which is simply a consequence of the dispersion relation and 
the larger inertial terms. The perpendicular group velocity 
remains mostly below 1 km/s, at least at the lower frequen- 
cies. 


5 Parallel electric fields 


the estimated perpendicular electric field ( 5E± ) and the prop- 
agation angle of the DAW through 

| «Jf?||(/) |= 9t i 5E x (f) | (8) 

The propagation angle of the DAW can be found from tan 6 = 
fc||/fcx. and is usually very small (0 < e = yfrnjmi). 
When a DAW propagates at larger angles, 8 > e, it is known 
as a slow ion-acoustic (SIA) wave (Seyler and Wahlund, 1996; 
Seyler et al., 1998; Wahlund et al., 1998, Stasiewicz et al., 
2000), The angle of propagation (Figure 4, panel d) is small 
for low frequencies. Above the local gyrofrequency for 0+ 
(near 30 Hz on orbit IV) the propagation angle reaches e, 
consistent with the proposition that SIA emissions start to 
dominate the ELF emission characteristics (Wahlund et al., 
1998). The 5E±/8B± ratio for a SIA wave should approach 
10 5 * * 8 — 10 9 m/s, which is also consistent with the observed 
values presented here (figure 3, panel g). 

The third method (Figure 4, panel c) of estimating &E\\ is 
based on equation (5) above. Parallel electric field strengths 
of up to 0. 1-1 mV/m are inferred from methods 2 and 3 (pan- 
els b and c) during active times, while significantly lower val- 
ues are obtained from method 1 . This discrepancy is ascribed 
to the difficulties in estimating a frequency dependent field- 
aligned current from the magnetic field measurements, i.e. 
method 1 is considered to be the most unreliable. Also, the 
increase at the upper frequencies in panel a is just the result 
of the fact that the fluxgate magnetometer reaches the noise 
level. In any case, the field-aligned electric field strengths are 
correlated rather well with the Poynting flux (Figure 3, panel 

f>. 

5. 1 Effect on DAW induced field-aligned current systems 


The parallel electric field associated with a DAW can be cal- 
culated in various ways. We shall adopt three methods and 
compare the results; one method makes use of the magnetic 
field measurements, the other two are calculated based on 
the perpendicular electric field and the above estimated wave 
vectors. 

The first method (Figure 4, panel a) is based on the esti- 
mated field-aligned current from the magnetic field measure- 
ments (6B±) along the satellite path ( Sx = v, c 5t) through 
the approximate relationship 


5 h = 


5B± _ 6B±_ 
IM)\± ^qSx 


( 6 ) 


The total field-aligned current fluctuations for orbit IV are 
displayed in Figure 5 (panel c). However, for the calculation 
at hand, <$?j| will be treated as frequency dependent through 
the magnetic field fluctuations. Subsequently the parallel 
electric field is calculated from the polarisation relation for 
inertial Alfv6n waves (DAW with cold electrons) as 


| 5E \\ (/) |= *)27r/A’ | | (7) 

The second method (Figure 4, panel b) is calculated from 


The total field-aligned current variation is calculated from 
the fluxgate magnetometer data according to equation 6 and 
displayed for orbit IV in Figure 5 (panel c). Values up to 200 
fiAlm 2 are reached. The largest current fluctuations, which 
are mainly due to the ELF electromagnetic emissions, are 
detected when the Langmuir probe (panel a) detects plasma 
density depletions/cavities together with enhanced Poynting 
fluxes (Figure 3, panel f). 

The net integrated current along the satellite path (Jj|) is 
displayed in panel b. This current gives information on where 
the large scale current systems are situated. The net inte- 
grated current over a cavity is small, as can be most clearly 
seen near +55 s, +90 s, +135 s, and +155 s; i.e., the effect 
of DAW activity is small on the overall field-aligned cur- 
rent system. Instead, the Alfv6nic activity creates a local 
region of field-aligned currents through the acceleration of 
the low energy plasma. Field-aligned electrons move down- 
wards by the parallel electric field and transverse ion heating 
occurs possibly through intermediate more electrostatic ELF 
waves. As a result the plasma density drops, which causes 
the increased perpendicular wave vectors and the inertial and 
kinetic dissipative effects enhancements. A small net field- 
alignedcurrent will remain due to the accelerated plasma. 
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Fig. 4. The field-aligned electric field as calculated by three different methods are presented in the three top panels (see text for further details), together with 
the estimated propagation angle of a dispersive Alfvgn wave (last panel). Parallel electric fields of the order of 0.1-1 mV/m arc inferred, which cm readily 
accelerate thermal electrons to energies of several 100 cV within fractions of a second. 


A larger integrated current can be found in the time in- 
terval 62-82 s (17:02:20 UT-17:02:50 UT), for which only 
low levels of DAW activity arc found, the downward pro- 
ton energy flux rises to 1 W/m 2 (Figure 1, panel b; Figure 
2 , panel a) and the precipitating electrons are characterised 
by an isotropic, smooth, inverted- V-like population (compare 
with Figure 1, panel c, and Figure 2, panel d). Two DAW 
active regions are situated at the edges of this region where 
the precipitating electrons probably give rise to the large in- 
tegrated current by the constant flux of isotropic electrons. 
Since the associated field-aligned electric fields in this mid- 
dle region arc probably very small according to commonly 
used currcnt-voltage relationships (e.g, Knight, 1973), this 
current, compared to that during strong Alfv 6 n activity, is not 
easily cancelled out by the cold dense background plasma. 


During periods when the ELF wave activity is enhanced 
without Langmuir probe detection of plasma density cavities, 
e.g. during tbe proton injections on orbits I, II , and III 
(not shown) and during the time period 120-130 seconds on 
orbit TV , there is also almost no net integrated current. This 
can be interpreted as follows. The Alfv 6 n waves accelerate 
particles due to the large inertial and kinetic components and 
by so doing try to create a small current However, since this 
occurs in a dense background plasma the current carried by 
the energetic electrons will quickly be cancelled out. 


6 Correlation study 

The electromagnetic ELF wave characteristics (wave ampli- 
tudes and SEx/SB ±) have been correlated with various elec- 
tron characteristics (energy, spectral shape, downward en- 
ergy flux) as well as the plasma density for five selected spec- 
tral wave emission frequencies between 1-46 Hz. The corre- 
lations are based on 16,600 points, i.e. the data presented in 
Figure 3. 

6 . 1 Particle spectral fits 

The electron spectral shape is approximated by an exponen- 
tial fit to the TESP data and characterised by the two param- 
eters 7 and 0 according to 

fo(E,<t>)='y(E) (9) 

where a represents the pitch angle and fo is the measured 
count rate. The electron spectral fits were carried out for 
pitch angles 0°, 45°and 90°for which a was set to 0, 0.1 and 
0.2, respectively, for simplicity. The parameter 0 represents 
tbe degree to which tbe electron population was isotropic, 
0 ~ 0 , or field-aligned, 0 « 0 (i.e. spectral shape), and pa- 
rameter 7 is the magnitude of the flux in the field-aligned di- 
rection. It is not possible to perform a lag correlation (shift in 
time) between the field fluctuations and various characteris- 
tics of the detected particles due to the fact that low frequen- 
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TIME (seconds) 


Fig. 5. The plasma density estimated from the Langmuir probe measurmaits (panel a) together with die estimated field-aligned current density assuming the 
Freja spacecraft encounters spatial structures (panel c), and the corresponding integrated net current (panel b). Alfvfri wave activity results in localised regions 
with rather small net integrated currents on top larger scale current systems. 


cies correspond to large times while the bursts occur during 
much shorter times. 

6.2 Panel set A: field- aligned count dependence 

The results of the parameter fits for an electron energy of 55 
eV are compared with the electromagnetic field parameters 
(Figure 6, panel set A). Contour lines of constant (tensity of 
points are drawn in the figure in order to resolve better the 
statistical shape of the correlation. 

Both the electric and magnetic field amplitudes, and con- 
sequently the 6E±/5B± -ratio, correlate positively with the 
electron flux counts (parameter 7), except for the higher fre- 
quency channels of the magnetic field fluctuations. From this 
it follows that the Poynting flux is also correlated with the 
electron flux (not shown). The decrease in correlation be- 
tween the field fluctuations at higher frequencies is due to 
the fact that the fluxgate magnetometer reaches its noise de- 
tection level, and the result is thus an artifact. 

6.3 Panel set B : electron spectral shape dependence 

Electron spectral shape values (|/3]> of less than 6 should 
be regarded as representing isotropic electron spectra, while 
\0\ > 6 indicate more field-aligned populations. The more 
negative the value; the more field-aligned the spectral shape, 
Figure 6. The more field-aligned electron populations appear 


when the electric and magnetic field amplitudes increase. 
The result is valid for all displayed frequencies. 

The electron spectra at the low energy of 55 eV show a bet- 
ter correlation with the amplitude of the waves compared to a 
higher electron energy (not shown). This is not so surprising 
since the low energy field-aligned electrons are suggested to 
be accelerated more easily by the Alfv6n waves than higher 
energy electrons (Andersson et al., 2000). Also, the high en- 
ergy electrons, which are observed with a broad pitch angle 
distribution (low value of /?), are suggested to arrive before 
the Alfv6n waves themselves. 

6.4 Panel set C: Electron energy flux dependence 

The total downward energy flux within the energy range be- 
tween ~ 20-50 eV is correlated with the field characteristics. 
The result from this correlation follows in much the same 
way as the correlation with the field-aligned electron count 
rate at 55 eV, but is better. The contour curves indicate a 
good correlation between both the electric and magnetic field 
amplitudes with increased downward energy flux. 

6.5 Panel set D; Plasma tensity dependence 

The correlations of the plasma tensity with the field ampli- 
tudes show that high plasma tensities (i.e. outside cavities) 
are not well correlated with the electromagnetic wave fields. 
Within plasma cavities, on the other hand, the electric field 
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Fig. 6. Correlations between various characteristics of the electrons and electromagnetic field fluctuations. Four types of correlations are presented; A-B 
correlate the field-aligned electron count rate at 55 eV and the electron spectral shape; C correlate the total downward electron energy flux below 50 eV; D 
correlates the plasma density with the field fluctuations. The field characteristics are represented by the electric, magnetic and 5Ex/5B± ratio at 5 selected 
frequencies (1, 2, 5, 14, and 46 Hz). 
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amplitudes become large. The magnetic field strength hardly 
correlates with the plasma density at all, indicating that the 
magnetic fluctuations are not important for the formation of 
plasma cavities and that it is the more electrostatic compo- 
nent of the ELF wave activity that is related to the cavity for- 
mation. Again, since the magnetic fluctuations are below the 
detection level of the fluxgate magnetometer at the highest 
frequency channel, the correlation loses its physical meaning 
at that frequency. 

The good correlation between wave activity and low en- 
ergy field-aligned electron precipitation, and the anti-correla- 
tion between wave activity and plasma density, has previ- 
ously also been reported by Chaston et al. (1999). 

7 Summary and conclusions 

In the four analysed events, wave-particle interactions oc- 
curred together with intense downward proton injection. The 
ion injection can therefore be assumed to be produced by a 
common process which also creates the Alfv6n wave activity, 
which in turn interacts with electrons through its dispersive 
induced field-aligned electric field in its wave front. Altema- 
tivley an ion injection gives rise to the Alfv6n wave activity 
which can accelerate electrons. A scenario can be visualised 
in which the ratio Ti/T e plays an important role by becom- 
ing large so that the Alfv6n wave approaches the electrostatic 
limit and begins to accelerate more and more electrons and at 
the same time transversely heats the bulk ion population and 
raises the Ti/T e ratio further (Clark and Seyler, 1999). This 
positive feedback process then increases the effect of the in- 
ertial and kinetic terms further, and more energy dissipation 
can occur and so on. 

The Freja mission was in operation during declining solar 
activity condition. The presence of Alfv£n waves is identified 
when the 5E±/5B± ratio becomes equal to or exceeds the 
local AIfv6n velocity. The variations of the plasma density, 
the Earth’s magnetic field strength and the ion composition 
with geocentric distance cause the Alfv6n speed to have a 
maximum below 2 R e geocentric distance. 

During orbit IV (6653), which is an extreme dayside event 
recorded by the Freja satellite, the largest plasma tempera- 
tures were encountered and density depletions/cavities were 
found. Also during this orbit, the 5E±/6B±-va\iQ became 
significantly larger than the local Alfv6n speed (ua) at fre- 
quencies above the local oxygen gyro frequency (about 30 
Hz). Furthermore, the largest perpendicular wave vector was 
found at these larger frequencies, and the wave character- 
istics became more electrostatic with increasing frequency. 
The data were therefore consistent with an interpretation in 
terms of dispersive Alfv6n waves Lysak ami Lotko (1996), 
which attain slow ion acoustic characteristics at larger fre- 
quencies (e.g. Seyler et al., 1995; Wahlund et al., 1998). 

The data presented in this paper support an earlier analy- 
sis (Chaston et al., 1999) that strong ELF wave activity, es- 
pecially in the electrostatic limit, is associated with larger 
scale plasma density depletions/cavities and at the same time 


enhances the low-energy electron fluxes along the geomag- 
netic field lines. The correlation study presented also gives a 
good correlation between electromagnetic field fluctuations 
and the presence of accelerated electrons with energies cor- 
responding to values below the local Alfv6n velocity. Intense 
electrostatic ELF missions were found to correlate well with 
the occurrence of density cavities. 

During periods of large plasma densities, broadband ELF 
wave activity and bursty electrons, the integrated field-aligned 
net currents were zero or very small. When larger scale 
density depletions/cavities were detected from the Langmuir 
probe measurements, the integrated field-aligned net currents 
became more significant although remained rather small. The 
fluctuating field-aligned current densities, however, are large, 
reaching several 100 /iV/m. The largest integrated field- 
aligned net currents were detected in association with peri- 
ods of low or zero ELF wave activity. At those times the 
fluctuations in the field-aligned current densities were on the 
contrary very low. A possible explanation of this behaviour is 
that the more isotropic electron energy distribution occurring 
during those times carries the current, not really affecting 
the local plasma, while the large parallel electric fields as- 
sociated with dispersive Alfv6n wave activity interact heav- 
ily with the local plasma, which in turn compensates for the 
fluctuating field-aligned currents induced by the impulsive 
Alfv6n waves. The flux tube associated with a large net cur- 
rent might be closed in the ionosphere, while the dispersive 
Alfv^n induced currents may have a more complex fate. 
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